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Aims. To better understand our Galaxy, we investigate the pertinency of describing the system of nearby disk stars in terms of a two-components 
Schwarzschild velocity distribution.Using the proper motion and parallax information of Hipparcos database, we determine the parameters 
characterizing the local stellar velocity field of a sample of 22000 disk stars. The sample we use is essentially the same as the one described by 
the criteria adopted to study the LSR and the stream motion of the nearby stellar population. 

Methods. The selected data is modeled with a two component Schwarzschild velocity distribution whose parameters were determined by a 
least-square regression. The celestial sphere was divided into 72 equal area regions used to determine the parameters minimizing the final 
velocity distribution function. 

Results. We verify that the results are not significantly different for the early type stars from the classical treatment using a single Gaussian 
population distribution. For late type stars in the subgiant branch, in contrast, we verify that the two-component model gives a much more 
satisfactory representation. Our results indicate for each spectral class the presence of these late type stars of a low velocity dispersion 
population, with cr„ ^ 20km.s"' coexisting with a high velocity dispersion population having cr„ ^ 40km. s"'. Both populations belong to the 
disk with scale heights of 120pc and 220pc, respectively, relative to the Galactic disk. 
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1. Introduction 

The kinematics of nearby disk stars is a longstanding issue in 
the literature has led to some major discoveries in the Galactic 
structure (see North 119901 for a complete review in this sub- 
ject). It is clearly important to properly characterize the motion 
of disk stars in the solar neighborhood. Most of our present 
knowledge of the dynamical properties, evolution, and origin 
of the Galactic disk depends on analyzing of the local stars. 
Hopefully this knowledge will also improve our understanding 
of the nearby galaxies where a similar star population is far 
more difficult to observe directly. A large body of information 
is already at our disposal to support our cuiTent knowledge of 
the disk (see for instance Binney & MeiTifield |1998l and Binney 
& Tremaine '1987' for a fully review of this subject). 

Normally the velocity field in the solar neighborhood is de- 
scribed by a single Schwarzschild (119071) velocity distribution, 
which is known to obey the collisionless Boltzmann equation. 
Due to the non-coUisional nature of the star fluid, the veloc- 
ity distribution can sustain different dispersions along the three 
Galactic coordinate axes. But this distribution does not remain 
strictly constant since the individual velocities of stars gradu- 
ally change on the Hubble time scale as a response to the dif- 
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ferent evolutionary processes involved. Therefore by studying 
the statistics of the presently observed velocities in the solar 
region, it is possible to infer the dispersions of the velocity dis- 
tribution and hopefully to see how they have changed with age, 
metallicity, and spectral class of the stellar population as pro- 
posed by many authors (Stromgren I 19871 Soubiran and Girard 
I2005l and references therein). 

A large body of data selected by age, metallicity, and spec- 
tral type attest to the validity of this bold description. One in- 
teresting issue that remains is how well we can separate the 
different kinematical populations of disk stars. It is a settled 
matter that disk stars of a given spectroscopic family share the 
same gross kinematical properties. Young early type stars are 
more concentrated towards the Galactic plane having a lower 
velocity dispersion in the range 10-20 km.s"'. In contrast late 
type stars, dominated by an older population, have a higher ve- 
locity dispersion (30-40 km.s ') and are more scattered around 
the Galactic plane (see Fig.|4]i. 

The same basic picture emerges when we sample nearby 
stars selected by their ages or metallicity. The young metal-rich 
population tends to have the lower velocity dispersion typical 
of the thin disk, while in the old metal-poor stars its value is 
definitely higher. These two large groups of stars also share 
the same velocity anisotropies, as measured by the ratio of the 
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velocity dispersions along the radial direction (ctu) and in the 
transversal direction ( both in the Galactic rotation direction, 
(Ty , and perpendicular to the plane of the disk, ctw ). 

A question that is also important in this context is whether 
an age spread should remain inside a given stellar population. 
Therefore we should expect that the young and the old stars of 
a given stellar type group should contribute differently to the 
overall dispersion velocities. Clearly for early type stars, the 
age spread is not expected to be large since these are short-lived 
objects, so we should also expect their kinematical properties 
to be more homogenous. But for late type stars we should be 
able to observe a much larger age spread, and therefore their 
global kinematics should be much more complex and sensi- 
tive to the age composition of the population. We should ex- 
pect for these objects a population of younger newly-born late 
type stars, therefore sharing the same kinematical properties as 
thin disk young stars, superposed on an old population of late 
type stars that have already suffered the cumulative effect of 
living inside a disk with gravitational potential irregularities. 
An alternative to this long-term, small gradual changing sce- 
nario was proposed by Carney et al. ("1989) based on possible 
detection of a discrete population of stars having [m/H] ^ -0.5 
and a velocity dispersion of 40km. s"' perpendicular to the disk. 
The existence of such a discrete population argues in favor of a 
relic population resulting from an early merger event occurring 
shortly after the formation of the disk. 

An earlier attempt to detect the effects of such a super- 
position of two different kinematical populations has already 
been done by Oort (I1932I I. who used a relatively small sam- 
ple of objects to conclude that this might be the case for late 
K-M dwarfs. According to that author, the kinematics of the 
two populations have cr„, - 10km. s"' and cr„. - 25km. s"'. In 
this paper we focus our attention on detecting a second kine- 
matical population using the proper-motion and parallax data 
of a sample of stars in the Hipparcos catalogue (ESA 1997). 
That is the best database available to date for that kind of in- 
vestigation because it provides positions, parallaxes, and an- 
nual proper motions for about 120 000 stars with an astromet- 
ric precision of Imas, allowing a discussion of this subject to a 
level of accuracy never attained before. Although the available 
radial velocity data is not complete, we can still use the homo- 
geneous set of proper-motion and parallax observations to de- 
duce the Galactic latitude and longitude components of the ve- 
locity field and constrain the kinematical description with that 
information. In other words, even if we cannot determine the 
individual components (U, V, W) for each star, we can develop 
a formalism to statistically determine the velocity dispersions 
CTj/, cry, and ctw for the whole population using this large sam- 
ple of objects. 

Among many others authors, Mignard (2000) has used this 
same approach to study the local kinematics. He constructed 
the Schwarchild probability distribution based on a selected set 
of approximately 22 000 Hipparcos stars by deducing the pa- 
rameters that characterize their local kinematic properties di- 
vided by spectroscopic groups. In the present work we use 
the same sample as Mignard(2000) to discuss another possi- 
bility for representing the probability distribution of velocities. 
Instead of the single Gaussian normally used to represent the 



Schwarzchild distribution, we used a two-Gaussian expression 
representing a two-fluid population. We used statistical tests to 
verify which of these descriptions gives a better representation 
of the real data. We present in Sect. 2 a discussion of the se- 
lection criteria to define a homogeneous sample of the nearby 
stellar population using the Hipparcos database. In Sect. 3 we 
present the basic description adopted to the local velocity field 
showing how the parameters of this distribution changes with 
the Galactic coordinate orientation. In Sect. 4 we discuss the 
kinematical parameters obtained in this work. Section 5 con- 
tains our main conclusions and prospects for future work. 

2. Material 

In Mignard (2000) we can find a very judicious selection of 
Hipparcos stars suitable for studying and discussing the local 
kinematics. He based his selection on five criteria that assert 
a sample of Hipparcos stars well-distributed in the celestial 
sphere and in spectral classes in a way that avoids any statistical 
biases (see Table[T]and Fig.[3]l. 

Table 1. Distribution of selected stars by spectral type interval 
estimated from the color index. 



B-V 


spec, type 


N,,„ 


0.00-0.15 


A0-A5 


4202 


0.15-0.30 


A5-F0 


3185 


0.30-0.45 


F0-F5 


2837 


0.85-1.15 


K0-K5 


6533 


1.15-1.40 


K5-M0 


3350 


1.40-1.60 


M0-M5 


2285 



The Mignard criteria, which we also adopted for selecting 
our sample from the Hipparcos catalogue, are summarized as: 
a) only single stars, due to their better astrometric mean pre- 
cision and better definition of their spectral type; b) stars lying 
between 0. 1 to 2.0 Kpc. The closer objects were eliminated due 
to their strong influence on the solution of the solar motion. In 
this selected range of distance, some of the individual paral- 
laxes might be not very precise, but we notice that even in this 
case their statistical meaning remains reliable for this study; 
c) through a detailed study of the completeness of Hipparcos 
stars, based on the more dense Tycho star catalogue (ESA 
1997), a completeness limit larger than 70% was established 
by a criteria using Galactic latitude, magnitude, and color; d) 
the stars were grouped in gross spectral types based on their 
colors; and e) stars with velocities higher than 60 - 90km.s 
depending on the spectral class, were excluded from the sta- 
tistical analysis. We retain the first four criteria in our sample, 
but for the last one adopt a less restrictive criterion discarding 
only those stars with peculiar velocity higher than 100. s ' . This 
would correspond approximately to a three sigma level cut us- 
ing the higher limit expectation for the disk velocity dispersion 
found in the literature (^ 30km. s '), so we feel confident that 
the high-velocity wing was faithfully sampled in our study. 

In Fig. [T] we present the latitude distribution of the sample 
selected by our criteria. The sharp contour of the selected ob- 
jects corresponds to the 70% completeness limit derived from 
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Fig. 1. Galactic latitude distribution of selected objects (blue points) and the rejected Hipparcos stars (green points). We can see 
the sharp contour corresponding to the 70% completeness hmit in relation to the Tycho catalogue. 



a comparison with the Tycho catalog. In Table[T]we present the 
distribution of the selected stars as a function of their spectral 
group. The sample is essentially identical to the one used by 
Mignard except for minor differences due to the fact that we 
adopted an analytical criterion describing the limiting magni- 
tude contour ( Vum = 7.68 + 1.15|(sin(/7)| - 0.43[(B - V) - 0.72] 
) to mimic the selection criteria based on the completeness. In 
particular, the color term was introduced to reproduce the vari- 
ation in the completeness limit with spectral type. The effect 
of the selection criteria on the HR diagram of our sample can 
be examined in Fig. |2] We can see that the early type objects 
in our sample are basically populated by young-disk main se- 
quence objects, and we exclude the very young ones that have 
peculiar kinematics associated to the region of their birth. In 
the case of late type objects the adopted criteria basically select 
subgiants and giants. All the lower main sequence objects were 
excluded by adopting the limiting distance and the magnitude 
detection limit of the satellite. The intermediate color objects 
(0.45 < B - V < 0.85) were also excluded due to the difficulty 
of assessing an spectroscopic type based only on the color in- 
formation. In this region, corresponding to the F5-K0 interval, 
the population is mixed by the presence of both dwarfs and gi- 
ant stars. 



The distribution of the selected objects in the celestial 
sphere is presented in Fig. [3] as expected, showing that early 
type stars are clearly more concentrated toward the plane of 
the Galactic disk. In contrast, late type giant stars are more 
uniformly distributed due to their larger scale height. This ef- 
fect is more clearly seen in Fig. H] where we show the den- 
sity of the distribution as a function of the vertical distance to 
the galactic plane. To obtain this density profile, we cut the 
selected sample into intervals of 100 pc of height above the 
Galactic plane counting all objects inside each of those ele- 
mentary cylinders. Early type stars have been combined into a 
single class A0-F5 and are tightly distributed around the plane 
of the disk. Late type objects have been combined into a single 
K0-M5 class and are clearly distributed over larger distances. 
For a single isothermal population of an autogravitating disk, 
we should expect that this density profile follow the relation 
Pdiz) = prfosech\z/zo) - p^o exp(-z2/2z2)^ where zo is the 
scale length and p^o represents the mass density in plane of 
the disk (Spitzer[T942]l. At moderate distances from the plane, 
this profile is approximately described by a Gaussian distribu- 
tion. The two continuous lines presented in Fig. |4] follow this 
density law and do not intend to really be a fit to the observed 
density distribution, but simply to reflect the eff'ect of the kine- 
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Fig. 2. HR diagram of the selected (blue points) and all Hipparcos stars (green points). 



matic properties determined in this work (see more details in 
Sects. 3 and 4) and serve to illustrate the scale length of these 
two populations. 

We can see that early type objects are relatively close to the 
thin disk, and their dispersion around the Galactic disk deter- 
mined directly from the density profile is cr. - 124pc. The den- 
sity profile of late type objects, on the contrary, is clearly spread 
out over larger distances as indicated by the height dispersion 
of cr. = 210pc. Nevertheless, both populations are clearly con- 
centrated inside the so-called Galactic thin disk as estimated 
by other authors (|z| < 350pc) (see for example Soubiran et al. 
12003 ). We note that the thick disk component would be visible 
only at larger distances from the Galactic plane (|z| ^ Ikpc) 
and is not sampled in our data, although there might be some 
degree of contamination. The point of maximum density on 
both curves presented in Fig.|4]does not correspond exactly to 
z = 0, meaning that the Sun is slightly displaced relative to the 
Galactic plane. 

3. Formalism 

The velocity vector of any star referred to the galactic center 
can be described as a composition of the local galactic rotation 



(V(R)) and a peculiar velocity (v,) related to the individual or- 
bit of that star in the galactic gravitational potential. The partic- 
ular value of these two vector fields depends on the stellar pop- 
ulation we are considering. Since these two components also 
affect the solar motion definition, it follows that the observed 
motion of any star relative to the Sun (v,) is also dependent on 
the stellar group considered 

V, = [V(7?)-y(7?o)] + [v.-vo]. (1) 

We refer the reader to Mignard (2000) for the use of this 
equation in connection to the Galactic rotation model obtain- 
ing the star's peculiar velocity. For a given homogeneous group 
of stars, the Oort constants were obtained by this author, and a 
careful analysis led to determining the local peculiar velocity 
field. In the present work we do not intend to repeat this anal- 
ysis, so we have opted to use his determination to recover the 
peculiar velocity field. Therefore by using the parallaxes and 
proper motions of the sampled stars, we could obtain the mo- 
tion of the Sun and the peculiar velocities necessary to the kine- 
matical analysis of the nearby disk stars. The components of the 
solar motion used here were extracted from Table 2 in Mignard 
(2000) for each interval of spectral type evaluated from the B-V 
color, as well as the Oort constants. 
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Fig. 3. Distribution of each spectral type in the celestial sphere. Early type stars are clearly more concentrated towards the Galactic 
plane while the late types are more evenly distributed. 



The peculiar velocity vector of a star in the Galaxy is deter- 
mined by the components u, v, and w, respectively directed 
toward the Galactic center (u), the direction of positive in- 
crease of Galactic longitudes (v), and the direction of the North 
Galactic Pole (w). For the purposes of our present analysis, in- 
stead of using this general dynamic system, it is rather more 
convenient to adopt a reference frame projecting these com- 
ponents along the hne of sight direction, Vr, and the two axes 
of Galactic coordinates, Vi and Vb. Since we are dealing with 
nearby objects, there is no loss of information, and the conver- 
sion between these two right-hand coordinate systems is easily 
made by adopting the transformations 



u = Vf cos / cos b -Vb cos / sin - vi sin / 

v= Vrsin/cosfo -Vbsin/sinfo -I- vicos/ (2) 



w = Vr sin b -i-Vb cos b. 

These equations enable us to express the nearby Galactic ve- 
locity field in terms of more directly observable velocity com- 
ponents that could be inferred from the parallax and proper mo- 
tion data. 

In general terms, the statistical distribution describing the 
Galactic structure of a group of stars in a specified unit vol- 
ume is defined in such a way that the number of objects in the 
velocity space dudvdw is 

6N = nfiu, V, w)d\idvdw = nf'{Vr, Vb< yi)dvrdvbdvi, (3) 

where the symbol n represents the local spatial density of stars, 
and the second identity is a direct consequence of the Jacobian 
of the above velocity transformations being equal to unity. In 
the present analysis, we need to deal with another Umitation: 




-800 -600 -400 -200 200 400 600 800 

z(pc) 



Fig. 4. Density of early and late type stars as a function of the distance to the galactic plane. The two curves represent the expected 
one component isothermal disk density distribution (p - pdosech^iz/zo)) with zo - 145.3 pc for early type stars (red continous 
line). For the late type star, the blue dashed line represents the expected composition of the two population having zo - 145. 3pc 
and Zo - 396.3 as explained in Sect. 4. 



the radial velocities are not known, at least for the vast major- 
ity of stars in our sample. Therefore we are restricted to using 
the velocity field as deduced from the parallax and proper mo- 
tion data, projected into the celestial sphere. For that reason it 
will be more convenient to work with the projected velocity 
distribution defined by 



6N = ng(wb, vi)dwbd\i. 



where the function 



gi^b, V/) 



X+Oi 
oo 



/'(V,., Vi, \i)d\r 



(4) 



(5) 



represents the velocity distribution integrated along the line 
of sight of the observer In writing this equation we implic- 
itly assume that our sample faithfully represents the radial ve- 
locity distribution. In that context we can define our statisti- 
cal problem as the determination of an expression for gi^b, v/) 
that could be applied to describing the Hipparcos parallax and 
proper motion database. 

As usual we adopt the view that the observed ellipsoid of 
velocities distribution of disk stars is aligned with the directions 



u, V, w and is approximately described by the so-called one- 
component Scwarzschild distribution 



/(u, V, w) 



where 



1 



(87r3)l/VuCrvCrv 



■ exp(e) 



1 u- v^ w2 

2 o-t o-t crt 



(6) 



(7) 



To ease the notation of the expressions below, we define 
the anisotropy parameters Fv = cTuIct^ and F„ - a-^lcr^, de- 
scribing the dispersions of the velocity ellipsoid in terms of the 
velocity dispersion towards the Galactic center. Moreover, we 
write the distribution function expressing the velocities in terms 
of the radial velocity dispersion (v^ = Vr/o-y, v/ = v//cru, v/, = 
Vfe/cu)- By adopting these definitions, we conclude that the 
transformation of the Swarzschild ellipsoid to the projected ve- 
locity distribution system is simply given by the expression 

Q = ~ ~ ^'^bb^l + CrlVr^l + Crbyr^b + ClbVlVb, (8) 
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where the numerical coefficients are easily determined by a di- 
rect comparison of this equation with the result of adopting the 
transformations of Eq.|2]into Eq|7] 



(9) 



Crr - cos^ / cos^ b + sin^ / cos^ b + rl, sin^ b 
cii - sin^ / + Fy cos^ I 

Chb - cos^ / sin^ b + rl sin^ / sin^ ^ + cos^ b 
Cri - (1 - rj) sin /cos Zcos 

Crb - (cos^ / + sin^ Z - r^) sin cos b 
cih - -(1 - Fy) sinZcos /sin/j. 

Using these coefficients we can obtain the corresponding 
projected distribution of velocities by an integration along the 
radial direction. 



•J —oo 



1 , 

exp[--c„.v,. + (ci-ivi + Crf,Vi)Vr]t/v,. (10) 



We can simplify this relation by using the the tabulated integral 



f 



exp(-/?^x^ + qx)dx — exp(-^) 



to obtain the expression 



givu Vb) 



r r 



rial 2 (CrfV/ + CriVi)2 

exp[--Q/v, - -CifcV^ + cibViVb + ] 



and, collecting the coeficients corresponding to the same power 
of V/ and Vb, we obtain the equation 



rvFvi; r 1 2 1 



rA 2 .yn^^^^~^1>>'^i - -^^hh^l + qibyiyb] (H) 



where the numerical coefficients are given by 



qn = 



F^ cos^ b + F^ sin^ I sin^ b + F^F?., cos^ / sin^ b 
cos^ I cos^ h + sin^ I cos^ b + VL sin^ b 



qbb 



qib - 



F^ (cos2 1 + Tl sin^ /) 



cos^ I cos^ b + T}, sin^ / cos^ b + sin^ b 

(F^ - 1 )F2 sin / cos / sin b 
cos^ I cos^ b + YI sin^ / cos^ + Fi. sin^ b 



(12) 



Depending on the choice of the observing direction, we ob- 
tain some particularly interesting cases for this projected ve- 
locity distribution. For example, in the direction = 0, we 
are observing the Galactic plane and therefore v^ = w. In 
that case the parameters describing the velocity distribution are 
qii - Til (cos^ / + Fy sin^ I), qi,i, - F^ and qu, - 0. Therefore the 
component of the velocity distribution along the latitude direc- 
tion Vb is identical to the distribution of the velocity ellipsoid in 



the direction orthogonal to the Galactic plane (w). Standing in 
the Galactic plane but looking now towards I - 0, the longitude 
velocity component (v;) is equal to the velocity component in 
the rotation direction v, while we have v/ - -u by looking at 

I = 7T/2. 

Another case of interest occurs by looking towards b = n/2 
at the North Galactic Pole. In that region qu - sin^ / + F^ cos^ I, 
qbb - cos^ I + rl sin^ /, and qu, - (F^ - 1) sin / cos I. If we ap- 
proximate the pole by keeping Z = 0, we have Vb - -u, v/ = v. 
And if I = n/2 we have v^ = -v, v/ = -u and qu = 1, qbb = T^. 

In principle, the bivariate distribution function g{vi, v/,) ob- 
tained can be directly compared with the observations at each 
direction on the celestial sphere. Since in each region the con- 
tribution of the fundamental vectors u, v, w changes with the 
Galactic coordinates, we should be able to deduce their in- 
fluence in the distribution function by comparing the predic- 
tion of this distribution with a large database. But unfortu- 
nately the number of stars in our sample is not that great 
for justifying the use of a fine grid, and we instead prefer to 
work with the marginal distributions i(v/) = J g(vi, yb)dvb and 
f(vfc) = / givi, yb)dvi and compare these distributions in some 
selected sectors of the sky. Using Eq. [TT]for the projected dis- 
tribution, we can easily obtain 



1 



s(v,) = -==exp(--^), 
\27To-i 2o-f 



(13) 



showing that the marginal distribution in Galactic longitude is 
Gaussian with a dispersion 



cos^ I 



(.- 



+ sin l)(T^ 



(14) 



that continously changes over the celestial sphere. 
Analogously, the marginal velocity distribution along the 
Galactic latitude is 



1 



, : exp( 

^J2ncrb 2cr, 



f(Vi,) = 



where 



cr? = (cos^ I sin^ b + 



sin^ / sin^ b 

f2 



COS 



(15) 



(16) 



Therefore the two marginal velocity distributions are also 
Gaussian, but their dispersions vary with the direction of obser- 
vation. The extension of these expressions to a two-component 
fluid system is straightforward, resulting in 

5N = n[a^i(vi, v,) -H (1 - a)g2{vb, vdjdvbdvi, (17) 

and the corresponding marginal distributions are 

i(v/) = asiivi) H- (1 - a)s2{-vi) 

f(Vfc) = Q'fl(V/,)-H(l-Q')f2(Vi) (18) 

where the dispersions of each component are given by Eqs.[T4l 
and[T6] and the parameter or is a measure of the relative impor- 
tance of the contribution of each component. 
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4. Analysis 

Since the dispersions involved in the marginal distributions de- 
pend upon combinations of square sine and cosine of / and 
b, we chose to compare these distributions with the observa- 
tional data by dividing the celestial sphere in 72 equal-area re- 
gions reduced to nine regions in the first octant of the sphere 
(0 < / < 90, Q < b < 90). For practical purposes, we adopted 
the three longitude segments located at 0-30°, 30°-60'', and 
60°-90''. The corresponding latitude segments were defined at 
b ^0,b ^ arcsin(l/3) ^ 19''.5, and b = arcsin(2/3) ^ 41°.8. 
By adopting that coarse division, each spherical sector has the 
same area, keeping a reasonable number of objects of 200-300 
in each region. Inside each region we binned the v;, v^, data in 
50 velocity intervals having Av = 4km. s"', therefore spanning 
the range of -100, -i-lOOkm.s ', which is adequate for studying 
the kinematics of the disk stars. 

A minimization code was projected to evaluate the 
sum of the square deviation of the marginal distributions 
s(cr/, cTy, Fy, r„; /, r(cri, cr„, F,,, F,,.; Z, /?). In each sector we 
consider that the dispersions cr; - cr/(/, b) and cr^ - cTbiJ, b) 
are constant, where I, b were evaluated at the centers of the cor- 
responding spherical sectors. Then, for each spectral type, the 
routine minimizes the combined mean square deviation of all 
the V/ and v;, histograms in the nine pre-defined celestial re- 
gions. Therefore the total number of bins is = 2x50x9 = 
900, and this database was used to estimate the mean square 
deviation. 

In Table|2]we present a summary of the determination of all 
the relevant kinematical parameters by the two ways: single and 
double Gaussian velocity distributions. For each spectral class, 
we have three rows of parameters. The first two rows indicate 
the minimization result using the two-component population 
model. The third row contains a determination of the param- 
eters based on a single kinematical population. In column A^, 
we indicate the total number of stars in each sample of spec- 
tral types, while < 5N^ >'^^ indicates the corresponding mean 
square deviation and its respective error estimated at a two- 
sigma level. In the next column we present the final chi-square 
per degree of freedom (x^/df). To estimate the chi-square, we 
evaluated the error at each velocity bin containing n objects as 
its Poisson expectation -\/n. In the other columns we have the 
proportion a in each case, the velocity dispersion in the u di- 
rection and the anisotropy parameters. For each parameter we 
also give the error bar evaluated at a two-sigma level. 

From the inspection of this table we can observe that, for 
early type stars (A0-F5), there is no significant net gain in de- 
scribing the kinematical data with the two-population model. 
For these objects the variation in the mean square deviation be- 
tween the two models is close to the expected intrinsic statis- 
tical fluctuation of their mean deviation. The same conclusion 
is reached when we compare the x^/df of the one- and two- 
component models. This one-component disk population can 
be directly compared with the results of Bienayme ( 19991 1 who 
obtained an average value of cr„ ^ 21. 3km. s"' for the early 
stars group, cr,, =^ 11. 3km. s"' (Fv, = 1.9),cr„, - 9.0km.s"' (F^, = 
2.4). Similar values were also obtained by Mignard (I2000I I who 



have determined cr„ = 19. 8km. s ',cr,, - 13. 3km. s '(F,, = 1.5), 
0-,, = 8.4km.s-'(F„ = 2.4). 

In Figs. |5] |6] and|2]we present the distribution of peculiar 
velocities for this group of early type stars in graphical terms. 
The top diagrams in these plots represent the results of the 
one-population model, while the lower panels shows the corre- 
sponding two-population model. The left panel represents the 
longitude velocity distribution, while the right panel represents 
the latitude velocity distribution. In both sets of panels, there is 
a small plot of the residual best model fit. In order to simplify 
the visual inspection we present these plots corresponding to 
the result of adding all the distributions for each spherical sec- 
tor in a single one. But we point out that the minimization code 
takes all the individually selected sectors into account since, as 
we saw above, the velocity dispersion varies with their position 
in the celestial sphere. In all these three groups correspond- 
ing to the early type stars, we can observe that the gain from 
considering a two-fluid population is barely marginal. In par- 
ticular the residual distribution has basically the same structure 
in all these cases, implying that there is no gain in considering 
a two-population model description. A more quantitative esti- 
mate was done using the Wald-Wolfowitz test to verify if the 
number of sign changes in the residue distribution is compati- 
ble with the null hypothesis of them being normally distributed. 
For the early type stars, there is no case where we could reject 
simultaneously both the v; and Vb distributions to the 95% con- 
fidence level. 

The more uncertain case among early type stars occurs 
in the A0-A5 class where the improvement in the fit (^(< 
6N^ >i/2) ~ 0.14) is barely larger than twice the error bar 
of the single-component solution (e = 0.05). Using the Wald- 
Wolfowitz test, we could reject the hypothesis that the residual 
V/ distribution is randomly distributed at a 95% confidence level 
but not V/,. Moreover the oscillatory behavior seen in Fig.|5]is 
removed when a low velocity dispersion population is added, 
giving some support to its reality. The effect is probable real, 
but we need a larger sample to be more positive about its de- 
tection. 

For the late type stars, our single population model basi- 
cally, as expected, recovers the same distribution of late type 
stars as found by Bienayme (I1999I I: cr„ ^ 34.4km. s"', cr,, ^ 
23.8km.s-' (F„ ^ 1.4), cr„ IT.Tkm.s"' (F^, ^ 1.9). And 
similar values were also found by Mignard (I2000I I. who esti- 
mated cr„ = 31. 2km. s"', cr,, - 22.0km.s"'(F,, - 1.4), cr„, - 
17.2km.s"'(F„, = 1.8). But contrary to what we found in early 
type stars, in the present case of the late type giant stars, there 
is an obvious improvement in the mean square deviation (< 
6N^ >i/2) by considering the two-fluid model description. For 
each spectral class, the improvement is greater than three times 
the mean expected error fluctuation, as we can verify in Table 
|2] In the one-population model we use the Wald-Wolfowitz test 
to reject the null hypothesis that the residual velocities are ran- 
domly distributed at the 95% confidence level in both v/ and 
Vh distributions in all the three spectral classes. The same type 
of improvement can be seen by the x^/df, which is systematic 
better in the two-fluid model. Therefore we feel that there is 
a real gain in describing this population by a two-component 
model. In gross terms the kinematics of late type giants can 




Fig. 5. Distribution of the longitude, left panels, and latitude, right panels, components of velocity. The two upper panels show 
the result of fitting a single population of the Scwarzschild velocity distribution. In the lower panels we present the result of 
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Fig. 6. Same as Fig.|5]for the A5F0 population. 
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Fig. 7. Same as Fig.|5]for the F0F5 population 
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Table 2. Kinematical parameters 



Sample N <SN^T^ YIW 



A0A5 4202 1.22±0.04 0.74 0.78±0.06 



19.54±1.60 
5.88±1.26 



1.63±0.33 
0.69±0.29 



2.74±0.30 
2.45±1.18 



1.36±0.05 



0.80 



15.25±1.07 1.35±0.21 2.65±0.30 



A5F0 3185 0.95±0.03 0.57 0.53±0.10 



24.85±2.88 
15.18±1.60 



1.64±0.46 
1.37±0.31 



2.58+0.59 
3.48±0.86 



0.99±0.03 0.57 



19.45±1.22 1.49±0.20 2.99±0.38 



F0F5 



2837 0.94±0.03 0.56 0.77±0.20 



24.50±2.25 
19.37±4.57 



1.59+0.28 
1.72+0.84 



2.19+0.43 
7.47+2.20 



0.96+0.03 0.55 



23.07+1.62 1.59+0.22 2.66+0.44 



K0K5 6533 1.01+0.03 



1.18 0.57+0.07 



40.66+4.48 
23.11+2.31 



1.49+0.37 
1.70+0.38 



2.07+0.51 
4.00+1.12 



1.19+0.04 



1.49 



31.37+1.60 1.60+0.22 2.84+0.47 



K5M0 3350 0.99+0.03 0.99 0.65+0.09 



41.25+5.19 
20.02+3.33 



1.54+0.47 
1.23+0.41 



2.07+0.59 
3.80+1.97 



1.09+0.04 



1.12 



31.57+2.84 1.42+0.26 2.56+0.54 



M0M5 2285 0.94+0.03 0.87 0.74+0.08 



41.28+5.20 
18.54+4.34 



1.32+0.36 
2.04+1.42 



2.37+0.70 
3.67+2.28 



1.06+0.04 



1.04 



— 33.10+3.50 1.43+0.33 2.70+0.65 



be described by a superposition of two well defined popula- 
tions. The low velocity dispersion component has an average 
dispersion ellipsoid (cr„,r,i, r„, = 21.1 + 3.3, 1.7+ .7, 3.8 + 1.2) 
that is similar to the average single-population solution found 
in early type stars (cr„,r,,r„, = 19.3 ± 1.3, 1.5 + 0.2, 2.8 ±0.4). 
Therefore it is tempting to interpret that these low velocity dis- 
persion populations as having the same origin as the thin disk 
population of early type stars. 

The high-velocity component of late type stars has a dis- 
tinct average ellipsoid:(o-„,r,.,r,, = 41.0 ±4.9, 1. 5 ±.0.4, 2.2 ± 
0.6) when compared to the early type population. The gain 
in the statistical description can also be clearly appreciated in 
Figs. [8] |9] and [TO] We observe from these plots that the sin- 
gle model population has the same distinct oscillatory resid- 
ual structure corresponding to a depression in the high-velocity 
wings of the distribution and an excess in the low-velocity cen- 
ter This occurs because the use of one single population is 
clearly not able to simultaneously reproduce the wings as well 
as the central low-velocity regime of the distribution. In the 
high-velocity wings, this one population model slightly under- 
estimates the observed data points, while the contrary occurs in 
the central region. 

The final result is that the minimum residual solution in this 
model results from a compromise trying to describe the low- 
velocity and high-velocity regions at the same time. This result 
is the oscillatory effect seen in the residual plots of the upper 
panels in these diagrams. The poorest results are obtained in 
bins where the residual reaches values in excess of 20 objects. 
Clearly the algorithm tries to find a compromise but is unable 
to deliver a randomly distributed error structure. In contrast, the 
consideration of a two-fluid population seen in the lower panels 



removes this oscillatory behavior, since the code uses one pop- 
ulation to describe the high-velocity wings and another one to 
describe the central low-velocity core. That is therefore why the 
final adjustment seen in Table|2]gives a much more satisfactory 
fit in the wings and also in the central region of low-velocity 
stars. The mean proportion of high-velocity dispersion stars, 
measured by the a parameter, coiTesponds to approximately 
65% of the objects, while the remaining 35% is composed of 
low dispersion stars. It is quite interesting that this rough pro- 
portion is maintained in all three spectral types with no clear 
trend toward systematic variation. 

In crude terms these two kinematical population found in 
the late type stars should have distinct scale heights above 
the Galactic plane in order to sustain their hydrostatic equi- 
librium in the gravitational field. In a multicomponent disk, the 
scale length of a given component is approximately given by 
Zo - cr^/^GI,od, where 'L„d represents the total projected mass 
density on the disk (Bahcall & Casertano. |I984| l. We note that 
the total mass density in the Galactic disk includes not only the 
stelar component but also the unseen contribution of the dark 
matter halo and other baryonic components such as gas and 
dust. 

The exact share of these components in the whole mixture 
is a matter of debate in the literature. According to Creze et 
al (I1998I I. the local dynamical mass deduced from Hipparcos 
data is po = 0.076 ± O.O15M0/pc-', and the total projected mass 
density should be close to Sq^/ - 33Mq/pc^, a lower value than 
previous estimates. Therefore, adopting this estimation, the two 
kinematical population found above should have a logarithmic 
scale length of zo = 54 pc and zo = 403 pc. The exact figure 
of these estimations is naturally subject to confirmation of the 
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Fig. 9. Same as Fig.|5]for the K5M0 population. 
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local disk density of matter. In Fig.|4]we have indicated the pre- 
diction corresponding to this estimative by the two continuous 
lines. In the case of early type stars, the estimation based on this 
model fits the observed density distribution observed from our 
sample quite nicely. For the late type giants, the observed dis- 
tribution deviates from this prediction probably because these 
stars are feeling the more extended disk potential that cannot 
be adequately be described by the thin disk approximation. But 
in both cases, it is clear that the two kinematical populations 
of late type stars can be considered as members of the Galactic 
thin disk. 

When we consider the whole range of spectral types, an- 
other interesting point is that the mean velocity dispersion, as 
deduced from the single population model, steadily increases 
by almost a factor of two when we consider the progression 
from the early type to the late type stars. As can be seen from 
Table |2l this progression of the velocity dispersion is quite 
smooth as we go from A0-A5 towards the M0-M5 class. But 
we can also see from our two-fluid model that the internal ve- 
locity dispersions of each late type population remains almost 
constant. Therefore the reason for the continous increase in the 
velocity dispersion is that the contribution of these two popu- 
lations changes as we go from the K0-K5 to the M0-M5 popu- 
lations. In the former case 57% of the stars belong to the high- 
velocity population, while almost 74% are high-velocity stars 
in the later case. Thus, even if the velocity dispersion of the 
populations of low-velocity and high-velocity stars remains ap- 
proximately constant, their combined dispersions increase con- 
tinously depending on their proportion. 

In fact, the two populations described above have quite 
similar properties to the dispersions of low and high metal- 
icity stars studied by Stromgren ( |1987l l. From this work the 
author concludes that the population of high-metallicity stars 
< [FejH] < -0.4 has a mean dispersion cr„ ^ 20km. s"', 
while the dispersion for low-metallicity -0.4 < [FejH} < 
-0.8 steadily increases reaching cr„ ^ 40km. s"' at the low- 
metallicity extreme. 

In Figs. [TT] and [121 we show the results for the two- 
population model applied to the combined class of the 
K0K5+K5M0+M0M5 sample. The v^ distribution in Fig. [H] 
shows a very good agreement of the model for all nine sec- 
tors. No systematic effect is present and the model represent 
all the sectors equally well. The v/ distribution shown in Fig. 
[TTI has a small systematic effect in the first three upper pan- 
els corresponding to the directions close to the galactic plane 
(b ^ 9. "8). The effect implies an excess of low-velocity stars 
in the database as compared with the prediction. In the wings 
of these three panels, we can also see a slightly oscillatory 
perturbation when the data points are compared to the model. 
Probably this effect is caused by a residual due to the subtrac- 
tion of the galactic rotational velocity that might be different for 
the two kinematical populations. The first upper panel of this 
figure should have a dispersion close to ctv, so at first glance, 
it would imply an excess of low-velocity stars in the direction 
of the Galactic rotation. But if that were the case, we should 
expect the same effect to be present in the right lower panel of 
Fig.[T2]and that is not the case. Therefore the possible origin of 
this excess is not clear to us. 



5. Conclusions 

Our analysis has shown that the kinematical properties of the 
population of early type stars (A0-F5) is described well by a 
single Swarzschild velocity distribution. No net gain is really 
detected when we use a two-population model even consider- 
ing the large number of objects in the present sample. Probably 
the detection of an additional minor population, if it exists, 
would require a much larger sample and a greater accuracy than 
the present one in the Hipparcos database, as expected by the 
Gaia experiment. There might be a hint of the presence of an- 
other small population since the residual plot on Figs. [TTI and 
[represents some structure in the low-velocity stars. But the 
size of the actual population was not large enough to permit a 
clear detection for our analysis. 

In the case of late type stars, the situation is clearly de- 
fined and points to the presence of two distinct kinematical 
populations. One population has high velocity dispersion stars 
with (Tu ^ 40km. s"', while the low velocity dispersion stars 
have cTy ^ 20km.s"'. The low-velocity component in the giant 
branch, corresponding to approximately 35% of this popula- 
tion, is typical of the thin disk kinematics. According to Caloi et 
al. ( |1999| see their Fig. 2), the youngest stars presently reaching 
the giant branch have ages in the range 5x10^ - IxlO^yr ,corre- 
sponding to objects with ages comparable to a few Galactic ro- 
tational periods. Therefore is tempting to identify these low ve- 
locity dispersion giants with thin disk stars that are now evolv- 
ing towards the giant branch. In that case the progression from 
the early type main sequence stars to these thin disk giants 
shows practically no variation in their velocity dispersion. In 
fact our result is consistent with a single population of thin disk 
objects having cr„, r,.,r„ - 20.2 + 2.3, 1.6 + .5, and 3.3 ± 1.1. 

The high-velocity late type giant population is typical of a 
more extended disk as discussed in several other analyses. It is 
quite improbable that our high-velocity sample could represent 
a transition to the thick disk. Gould et al. (1996) found that the 
solar neighborhood contains 80% of stars belonging to a disk 
having a scale length of 326 pc, while the other 20% of the 
population belongs to an extended disk with a scale length of 
656 pc. Since we are sampling a scale length well inside this 
thin disk, the degree of contamination by the more extended 
thick disk should be even lower than this figure. According to 
Soubiran ( |19931 l, this thin, old disk population is characterized 
by a-jj = 44 ± 6km. s"' and o-y - ± 4km. s"', close to the 
solution we found for our higher velocity giants. On the con- 
trary, the thick disk would require a higher velocity dispersion 
of o-jj =56+1 1km. s"' and crj/ = 43 + 6km. s"'. Although we 
cannot discard the presence of some degree of contamination, 
our feeling is thus that our high velocity dispersion late type 
population belongs to the thin disk. 

Another piece of information comes from the work of Caloi 
et al. (1 19991 1, where they study a large sample, including both 
giants and main sequence stars, using Hipparcos parallaxes and 
radial velocities extracted from the literature. An interesting 
point raised by their discussion is that stars having low space 
velocity relative to the LSR (-40km.s"' <V < lOkm.s"') tend 
to present a higher proportion of young objects (age< lO^yr). 
On the other hand, the bulk of objects in their sample having 




Fig, 10. Same as Fig.|5]for the MOMS population. 
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Fig. 11. Comparison of the longitude velocity distribution of the two models with the combined data of K0K5-I-K5M0-I-M0M5 
late type stars. The panels correspond to the nine sectors in the celestial sphere seen in Galactic coordinates. The three upper 
panels correspond to the the Galactic plane (t> ^ 9". 8), the second row are the intermediate zone (b - 30°.?), and the lower 
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Fig. 12. Same as Fig.[TT]for the latitude velocity component. 
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velocities comparable to the high-velocity tails of our sample 
(-60km.s ' < V < -40km.s ') appear older than lO^yr. By 
increasing the velocity range, they begin to sample the transi- 
tion zone to the thick disk and the inner halo objects. Therefore 
our present result supports the idea that our high-velocity pop- 
ulation is due to local disk stars having ages comparable with a 
few rotation periods of the Galaxy. In that case the dynamical 
differentiation detected in this work could probably reflect the 
diffusion of stellar orbits due to the dynamical relaxation due 
to the irregular time-varying gravitational field in the Galactic 
disk (Wielen [T977] l. In particular, Griv et al. (120011) estimate 
that scattering by potential irregularities in the galactic disk 
could cause a velocity variation in lO^yr comparable to the dif- 
ference we have detected here. Therefore, in this scenario, the 
older objects would correspond to the population of the high- 
velocity tail of our late giant distribution. 

There is however one point that might challenge this inter- 
pretation, since our high-velocity population present no sign of 
variation among the spectral classes K0-K5, K5-M0, and MO- 
MS. In all of these three samples the results are remarkably 
consistent with a single kinematical population that is quite dis- 
tinct from the one found in the low-velocity thin disk. There is 
no sign of a smooth transition between the two regimes, which 
would be expected in the presence of a gradual migration phe- 
nomenon such as the scattering by irregularities. In that case 
one could naively expect that the kinematical properties would 
smoothly change as we sample objects with different classes 
and probably different ages and chemical composition. That is 
not what we see since the kinematical properties of our sam- 
ple of late type giants are quite consistent with a single kine- 
matical population. In that context it is interesting to note the 
work of Norris ( 1 19871 1, according to which the prevalence of 
the scattering mechanism implies a smooth transition regime, 
whereas the presence of two overlapping distributions would 
favor a discrete-origin model. In fact, Carney et al. 119891 use 
this argument in favor of the existence of a discrete popula- 
tion in the disk with [m/H] ^ -0.5 and a velocity dispersion 
perpendicular to the disk of 40kms"'. It is suggested that this 
population could be the signature of an early merger event oc- 
curring shortly after the disk had formed. 
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